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Abstract. A method is developed for realizing entangled states and general second 
quantized fermionic and bosonic fields in the framework of the fermionic projector. 
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1. Introduction 

In [3l |9] it was proposed to formulate physics based on a new action principle in 
space-time. One difference of this approach to standard quantum field theory is that 
a many-particle state no longer corresponds to a vector in the fermionic Fock space, 
but instead it is described by the so-called fermionic projector, an operator which acts 
on the one-particle Hilbert space (or more generally on an indefinite inner product 
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space spanned by the one-particle wave functions). Another difference is that the 
bosonic fields obtained in the so-called continuum limit are only classical. Due to 
these differences, it is not at all obvious whether the fermionic projector can account 
for all quantum effects observed in nature. More specifically, is it possible to describe 
entanglement? Can one reproduce the effects of second quantized fields? 

In this paper, we shall analyze these questions in detail. We will show that it is in- 
deed possible to describe entangled states as well as general second quantized bosonic 
and fermionic fields in the framework of the fermionic projector. This is achieved by in- 
troducing the physical concept of a microscopic mixing of decoherent subsystems. The 
physical picture is that space-time is not smooth on the microscopic scale (typically 
thought of as the Planck scale), but has a non-trivial microstructure. Homogenizing 
this microstructure, we obtain an effective description of the system by a vector in 
the Fock space. To make this picture precise, we use the fact that in the framework 
of the fermionic projector, the usual topological and causal structure of Minkowski 
space is not a-priori given, but it is induced on the space-time points by the states of 
the fermionic projector. Thus by bringing the wave functions between certain pairs 
of space-time points "out of phase" , we obtain decoherence effects which result in a 
decomposition of the whole system into subsystems between which the usual causal 
relations are no longer valid. This makes it possible to realize many independent phys- 
ical systems simultaneously in one spacetime, in such a way that homogenizing on the 
microscopic scale leads to an effective "superposition" of the subsystems. For techni- 
cal simplicity, we will describe the microscopic mixing by localizing the subsystems in 
disjoint space-time regions (see Figure [2] on page I16p . But one can also think of the 
subsystems as being delocalized, similar as if one combines several images in a single 
hologram (see Section l5T3|) . 

In order to make the paper self-contained easily accessible, we begin in Chapter [2] 
with a brief outline of the fermionic projector approach. In Chapter O we recall the 
basics on entanglement and Fock spaces and work out the connection to projectors 
in the one-particle Hilbert space. In Chapter |4] we describe entanglement using the 
concept of microscopic mixing and by introducing a decoherence between space-time 
regions. As is worked out in Chapter [5j this notion of decoherence also makes it 
possible to describe second quantized bosonic fields. 



2. An Outline of the Fermionic Projector Approach 

In this chapter we give a brief introduction to the framework of the fermionic pro- 
jector, outlining a few ideas, methods and results. For details we refer to [31 El [8] or 
to the review papers [101 [9] . 

2.1. An Action Principle for Fermion Systems in Discrete Space-Time. In 

the fermionic projector approach, the physical equations are formulated intrinsically 
in a discrete space-time. To introduce the basic framework, we let (H, <.|.>) be 
a finite-dimensional complex vector space endowed with an indefinite inner product 
(thus <.|.> is non-degenerate, but not positive definite). Next, we let M = {1, . . . , m} 
be a finite set. To every point x G M we associate a projector E x (a projector in H 
is defined just as in Hilbert spaces as a linear operator which is idempotent and self- 
adjoint). We assume that these projectors are orthogonal and complete in the sense 



ENTANGLEMENT AND SECOND QUANTIZATION 



3 



that 

E x E y = 5 xy E x and ^ E x = 1 . (2.1) 

Furthermore, we assume that the images E X (H) C H of these projectors are non-dege- 
nerate subspaces of H, which all have the same signature (2, 2). The points x S M are 
called discrete space-time points, and the corresponding projectors E x are the space- 
time projectors. The structure (H, <.|.> , (E x ) xe M) is called discrete space-time. The 
particles of our system are described by one more projector P in H, the so-called 
fermionic projector, which has the property that its image P{H) is a negative definite 
subspace of H. The resulting system (H, <.|.> , (E x ) x ^m, P) is referred to as afermion 
system in discrete space-time. 

Let us briefly discuss these definitions and introduce a convenient notation. The 
vectors in the image of P have the interpretation as the occupied fermionic states of 
our system, and thus the rank of P gives the number of particles f := dim P(H). The 
space-time projectors E x can be used to project vectors of H to the subspace E X (H) C 
H. Using a more graphic notion, we also refer to this projection as the localization 
at the space-time point x. For the localization of a vector tp G "K we use the short 
notation 

%b(x):=E x ^ (2.2) 
and refer to ip{x) as the corresponding wave function. Having the connection to rela- 
tivistic quantum mechanics in mind (see Section [2.31 below), we refer to E X (H) as the 
spinor space at x and denote it by S x . It is endowed with the inner product <.\E X .> 
of signature (2,2), which we also denote by -<.|.>-. Using the relations (12. ip . we can 
then write 

<t/#> = ]T -<il>(x)\<f>(x)>- . (2.3) 

The localization of the fermionic projector is denoted by P(x,y) := E x PE y . This 
operator maps S y C H to S x , and we usually regard it as a mapping between these 
subspaces, 

P(x,y) = E x PE y : S y -> S x . 
Again using (|2,ip . we can write the wave function corresponding to Pip as follows, 

(P^)(x) = E X P^ = Y.ExPEyiP = J2( E * PE v)( E vTp)> 

j/eM j/eM 

and thus 

(PiP)(x) = £P(z,j/)^(y). (2.4) 

yeM 

This relation resembles the representation of an operator with an integral kernel, and 
therefore we call P{x, y) the discrete kernel of the fermionic projector. Finally, it is 
often useful to choose an orthonormal basis ipi, . . . ,ipf of P{H) (i.e. <ipi\ipj> = — <5«). 
Then the fermionic projector and its discrete kernel can be written in bra/ket-notation 
as 

/ / 
P = -J2^ >< ^\ and P(x,y) = -Y,\^( X >^M\ ■ ( 2 -5) 

In order to introduce our action principle, for any x, y £ M we define the closed 
chain A xy by 

A xy = P(x,y)P{y,x) : S x ^ S x . (2.6) 
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We denote the eigenvalues of A xy counted with algebraic multiplicities by \%, . . . , A4 
and define the spectral weight \A xy \ by 

4 

Similarly, one can take the spectral weight of powers of A xy . Our action principle is 
to 



minimize the action S[P] := \A xy \ 



(2.7) 

under the constraint T[P] : = |^-xy| 2 = const , 

x,yeM 

where we consider variations of the fermionic projector, keeping the number of particles 
as well as discrete space-time fixed. In [6] it is shown that minimizers of this nonlinear 
variational principle exist. For a discussion of the underlying physical principles see [U 
Section 2]. 



2.2. The Correspondence to Minkowski Space. At first sight, the above setting 
seems inappropriate for physical applications because important structures like the 
notion of causality, topology and metric of space-time, gauge fields, etc., are miss- 
ing. However, the idea is that these additional structures arise as a consequence of 
a self-organization of the states of the fermionic projector as described by our action 
principle. More specifically, for a given minimizing fermionic projector P, its dis- 
crete kernel P(x, y) induces relations between the points x and y, which can be used 
to introduce additional structures in space-time. This mechanism is referred to as 
spontaneous structure formation. The first rigorous result on spontaneous structure 
formation was obtained in [5], where it is shown that the permutation symmetry of 
the space-time points is spontaneously broken, giving rise to a non-trivial outer sym- 
metry group. More detailed information is obtained by analyzing the eigenvalues of 
the closed chain: 

Definition 2.1. (causal structure) Two space-time points x,y 6 M are called time- 
like separated if the spectrum of A xy is real. The points are spacelike separated if the 
spectrum of A xy forms two complex conjugate pairs having the same absolute value. In 
all other cases, the two points are lightlike separated. 

Furthermore, from the discrete kernel P(x, y) one can deduce objects of differen- 
tial geometry like the spin connection and curvature (see jllj). Finally, the pa- 
pers [7J [12] show that in the limit of an infinite number of particles and space-time 
points, our action principle admits minimizers which are regularizations of vacuum 
Dirac sea structures in Minkowski space (see also |10j). In this limit, the wave func- 
tion i/j(x), (12. 2|) . goes over to a four-component Dirac spinor. The spin scalar product 
becomes -<if:{x)\(j){x)>- = ip(x)4>(x), where ip = ■0'7° is the usual adjoint spinor of 
Dirac theory. Discrete space-time goes over to a space-time continuum M = K , and 
the sum in (|2.3p becomes a space-time integral, 



<ip\4>> : = / ^.i(j(x)\4>(x)>- d 4 x 
Jm 



(2.8) 
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In the simplest case of one sea, the discrete kernel corresponds to the Fourier integral 
of the lower mass shell, 

(% 7 j + m) 5(k 2 - m 2 ) Q(-k°) e ~ lk ^ , (2.9) 

where is the Heaviside function (more generally, one can take sums or direct sums 
of Dirac seas to describe different types of elementary particles; see §5.1]). Here 
the exponent e~ ik \ x ~y) involves the Minkowski metric. Even more, (|2.9p determines 
the Minkowski metric and can be used define it. If this is done, the causal structure 
corresponding to the Minkowski metric indeed agrees with the spectral definition in 
Def. 12.11 (for details see [H Section 3]). Moreover, one can introduce all the familiar 
objects of Dirac theory. For example, the non-negative quantity -^tp\ r y°ijj>- has the in- 
terpretation as the probability density of the Dirac particle. Polarizing and integrating 
over space yields the scalar product 

(V#)=/ ^(t,x)\j°0(t,x)ydx. (2.10) 

J i=const 

For solutions of the Dirac equation, current conservation implies that this scalar prod- 
uct is time independent. 

2.3. The Continuum Limit, a Formulation of Quantum Field Theory. The 

above correspondence to vacuum Dirac sea structures can also be used to analyze our 
action principle for interacting systems in the so-called continuum limit (for details 
see [U Chapter 4] and [8]). We now outline a few ideas and constructions needed 
later on. First, it is helpful to observe that the vacuum fermionic projector (|2.9p is a 
solution of the Dirac equation (ij-'dj — m)P sea (x, y) = 0. To introduce the interaction, 
we replace the free Dirac operator by a more general Dirac operator, for example 
involving gauge potentials or a gravitational field. Thus, considering for simplicity an 
electromagnetic potential A, we demand that 

(i-y j (dj - ieAj) - m) P(x,y) = . (2.11) 

Moreover, we introduce particles and anti-particles by occupying (suitably normalized) 
positive-energy states and removing states of the sea, 

-, n f , n a 

P(x, y) = P sca (x, y) - — J2 \M*><My)\ + y E \Mx>*Mv)\ ■ (2-12) 

k=l 1=1 

Using the so-called causal perturbation expansion and light-cone expansion, the fermio- 
nic projector can be introduced via (12. lip and (12.121) . 

It is important that our setting so far does not involve the field equations; in particu- 
lar, the electromagnetic potential in the Dirac equation (|2.1ip does not need to satisfy 
the Maxwell equations. Instead, the field equations should be derived from our action 
principle (|2.7p . Indeed, analyzing the corresponding Euler-Lagrange equations, one 
finds that they are satisfied only if the potentials in the Dirac equation satisfy certain 
constraints. Some of these constraints are partial differential equations involving the 
potentials as well as the wave functions of the particles and anti-particles in (j2. 12|) . 
In [8], such field equations are analyzed in detail for a system involving an axial field. 
In order to keep the setting as simple as possible, we here consider the analogous field 
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equation for the electromagnetic field 

n f n a 

d jk A k - UAj = eJ2^khji>k>~ ~ e^2^i\jj(/>ty . (2.13) 
fe=i l=i 

With (12. lip and (12.13p . the interaction as described by the action principle (12. 7p re- 
duces in the continuum limit to the coupled Dirac- Maxwell equations. The many- 
fermion state is again described by the fermionic projector, which is built up of one- 
particle wave functions. The electromagnetic field is merely a classical bosonic field. 

For the considerations in Chapters 0] and [5l it is important to keep in mind that in 
the framework of the fermionic projector, space-time is not smooth on the microscopic 
scale, but it has an underlying discrete structure. The dynamics is described intrin- 
sically in discrete space-time by the action principle (|2,7p . A minimizing fermionic 
projector has a rich microscopic structure from which one can deduce notions which 
have a correspondence to macroscopic physics. In the continuum limit, the causal and 
metric structure of Minkowski space can be recovered from the fermionic projector 
using the notion of Definition 12.11 Thus we can say that the wave functions of the 
Dirac particles (also taking into account the states of the Dirac sea) generate the causal 
and geometric structure of space-time. This observation will be helpful in Chapter [H 
where by bringing the wave functions in different regions of space-time "out of phase" , 
we will be able to turn off causal influences between these regions. 

3. Preliminaries on Projectors, Fock Spaces and Entanglement 

In this chapter we first recall the notions of the fermionic Fock space and entan- 
glement, also fixing our notation. Then we show that a projector in the one-particle 
Hilbert space uniquely determines a Hartree-Fock state in the fermionic Fock space, 
making it impossible to describe entangled states. 

3.1. The Fermionic Fock Space and Entanglement. In non-relativistic quantum 
mechanics, the one-particle states form a separable Hilbert space (IK, (.|.)). Simi- 
larly, in Dirac theory we let "K be the Hilbert space corresponding to the scalar prod- 
uct (|2,10p . In each of these settings, a many-fermion state is usually described by a 
vector in the fermionic Fock space, which we now introduce (see also |17} Section II. 4] 
or [TH Section 1.1]). We let "K n = % <g> • • • <g) Oi be the n-fold tensor product, endowed 
with the natural scalar product 

(i>i ® • • • ® i> n I <h. ® • • • ® <M ■= (V>il<M • • • {i>n\4>n) ■ (3.1) 
Totally anti-symmetrizing the tensor product gives the wedge product 

ih A • • • A := ^ (-l) Sign(<T) ® • • • ® VV(n) (3.2) 

(here S n denotes the set of all permutations and sign(cr) is the sign of the permuta- 
tion a). The wedge product gives rise to a mapping 

A n : IKx. . .xJK -»• !K n : (V>i, • • • , ip n ) ^ ipl A • • • A tp n . 

s v y 

n factors 

We denote the image of this mapping by 3^ F . The vectors in 3~ HF are called n- 
particle Hartree-Fock states or factorizable states. These states do in general not form 
a vector space, as the following example shows, which in discussions of spin correlation 
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experiments and Bell's inequalities is often referred to as the EPR singlet state (see 
for example [TJ Section 1.5]). 

Example 3.1. (The spatially separated singlet state) We consider the one-particle 
Hilbert space J£ = C 2 A C| of two spins (in quantum information theory called 
"qubits"), located at the positions of two observers Alice and Bob. Choosing in C 2 
the standard basis ^ = (1,0) and V ; = (0, 1) yields the basis (4>4>4>4) of <K - 
The spatially separated singlet state is the following linear combination of 2-particle 
Hartree-Fock states 

* := ^ (VWWWL) • (3-3) 

Let us verify in detail that this state is not factorizable. Thus assume conversely that 
can be written as a product, 

^ = ipi A ip 2 ■ 

Computing this wedge product in the basis representations 

ipi = 44 + 44 + "b^b + 44 (3-4) 

^ 2 = 44 + 44 + 44 + 44 , (3-5) 

the vanishing of the term ~ ya^A i m pli es that the vectors a A ifj A +a A ip A and P A ip A + 

ii *f i 

(3 A ip A must be linearly dependent. Similarly, the vanishing of the term ~ ip B A tp B 

implies that the vectors a B ip B + a B w B an d Pb^b + Pb^b are hnearly dependent. 
Hence ip2 can be written as 

ip 2 = Pa (44 + 44) + Pb (44 + 44) ( 3 -6) 



with suitable complex coefficients /3a and Pb- Taking the wedge product of (|3T 
and ()3.6p yields 



* = Vi A V2 = (Pb - Pa) (aVl + A (44 + 44 



Multiplying out and comparing with (13. 3D . one sees that the products 44 and 44 
must be non-zero, and thus none of these four coefficients vanishes. But then the 
term ~ A ip B is non-zero, a contradiction. 

We denote the vector space generated by the n-particle Hartree-Fock states by 



J n = <A n (^™)> . 
Their direct sum is the fermionic Fock space, 



n=0 



The non- factorizable vectors Vl/ 6 3~ n \A n (J{ n ) are called entangled states. The spatially 
separated singlet state is the standard example of an entangled state. Entanglement 
is a basic phenomenon of quantum physics with important potential applications in 
quantum computing. 
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3.2. Projectors and Hartree-Fock States. Let us examine in which sense a pro- 
jector in the one-particle Hilbert space characterizes a many-particle quantum state. 
Thus let P be a projector in the Hilbert space (JH, {■{■)), for simplicity of finite rank /, 
i.e. 

p* = p = p 2 and dimP(5£) = /. (3.7) 
In order to get a connection to the fermionic Fock space formalism, we choose an 
orthonormal basis fa, . . . , ipf of P(IK) and form the Hartree-Fock state 

* := fa a • • • A V/ G • (3.8) 
The choice of our orthonormal basis was unique only up to the unitary transformations 

/ 

fa^fa = J2 U ij^j with Ue\J{f). (3.9) 

3=1 

Due to the anti-symmetrization, this transformation changes the corresponding Hartree- 
Fock state only by a phase factor, 

fa A • • • A 4> f = detU fa A • • • Aip f . (3.10) 

Thus we can indeed associate to the projector P a Hartree-Fock state, which is well- 
defined up to a phase. As the phase of has no physical significance, the physical 
system is described equivalently by a projector Pf on the many-particle state , i.e. 
in bra-/ket notation^ 

P / = ^yp-|*)(*| = /!|^lA---AV/)(^iA---A^ / | :S f ^W f . (3.11) 

ll 11^ 

Since the phase freedom drops out when forming the projector fj3. 1 1 f) . this operator 
is well-defined. The next proposition gives an alternative definition of Pf which does 
not involve a choice of basis. 

Proposition 3.2. For any projector P in (IK, (.|.)) of rank f , the corresponding op- 
erator 

P f : 2t S -*2r s : fa A • • • A ip f -»■ (Pfa) A • • • A (Pip f ) (3-12) 

is a projector onto an f -particle Hartree-Fock state. The mapping P — )• Pf gives a one- 
to-one correspondence between projectors in VL and projectors on Hartree-Fock states 
in 3~. 

Proof. It follows immediately from the definitions that Pf is symmetric and idempo- 
tent, and is thus a projector. To compute the rank of Pf, we choose an orthonormal 
basis fa,...,ipf of P(!H) and extend it to an orthonormal basis of "K. As is obvi- 
ous from (|3.12p . the operator Pf applied to any wedge product of basis vectors van- 
ishes unless all basis vectors are elements of the set {fa, . . . ,ipf}. Hence the vector 
^ := fa A • • • A ipf is a basis of the image of Pf. We conclude that Pf has indeed rank 
one and is thus a projector onto the Hartree-Fock state ^. 

^Here the factor /! comes about because, according to our conventions (|3.2p and (|3.1|l . 

(V>i A---Aipf\ipiA---A i>f) = {V>i A • • • A ipf | Vi &> • • • ® 4>f) 

= ji E (-i) sisn((T) WwdIiM • • • Wv(/)IV7> = 4 • 
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Now suppose conversely that Pf is a projector onto a Hartree-Fock state. Repre- 
senting this operator in the form (|3.1ip . we let P be the projector in IK on the sub- 
space <ipi, ■ ■ ■ , ipf> ■ Then the operator Pf has the representation (|3.12p . concluding 
the proof. □ 

3.3. Projectors and Expectation Values. We now consider how expectation values 
of observables can be expressed in terms of the projectors P and Pf. We begin with 
a one-particle observable O, being a self-adjoint operator in the one-particle Hilbert 
space "K. By 

J (V>i A • • • A ip n ) := (Ctyl) A • • • A if; n + ip! A (Ofa) A • • • A ip n 

+ • • • +ipl A ■ • • A ipn-l A (Oip n ) 

we can define a corresponding operator O 3 " on the Fock space 3". This operator pre- 
serves the number of particles in the sense that it maps the n-particle subspace 3~ n to 
itself. 

Suppose that an /-fermion state is described by a projector P, (|3.7p . The next 
lemma shows how the expectation values of O 3 " and of products of one-particle oper- 
ators can be expressed in terms of traces involving the projector P. 

Lemma 3.3. Suppose that O and 0\n are one-particle observables. Describing a 
many-fermion state by a projector P in (JC, (.|.)), we have 

(O 37 ) = Tr M (PO) (3.14) 

(Of Of) = Tr^POxOa) + Tr^PO^ Tr K (P0 2 ) - Tr Ji (P0 1 P0 2 ) . (3.15) 

Proof. The expectation values are obtained by taking the trace of the observables 
multiplied by the operator Pf, 

(O j ) = Tr ?/ (P/O j ) , (Of Of) = Tr^PyOfOf) . 

Representing Pf in the form (13. IIP , it follows that 

(0 J ) = f\ (ipx A ■ ■ ■ Aifif \ O y (V>i A • • • A i> f )) 

= /! (tpx A • • • A ipf | O 3 (ip! ® • • • ® ipf)) 

f 

<t£S,i i=l 

where in the last step we applied (|3.13p together with (|3.ip and used the fact that the 
vectors ip\, . . . , ipf are orthonormal. This proves ()3.14p . 
With the same method, we obtain 

(Of Of ) = /! A • • • A ip f | Of Of (^i ® • • • ® V/)) 

/ 

= /! (V>i A • • • A V/ I V>i ® • • • i/> k -i ® (Oi0 2 ^ fe ) 8) V'fc+i • • • V'/) 
fc=i 

+ /! A • • • A ipf | Vi ® • • ■ (Oiipk) ■ ■ ■ (0 2 ipi) ■ ■ ■ ® ipf) 

k^l 

fc=l fc^Z 
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The last sum can be extended to all k, I = 1, ... , /, because the summands for k = I 
vanish. We thus obtain (|3. 15[) . □ 

The method of this lemma immediately extends to higher powers of one-particle ob- 
servables. 

More generally, one can consider many-particle observables, described by a self- 
adjoint operator O on the Fock space 3". In this paper, we shall only consider observ- 
ables which preserve the number of particles, i.e. which are invariant on the n-particle 
subspaces 3~ n , 

O : J n -> J n . (3.16) 

This assumption is justified by the physical law of the conservation of the baryon 
and lepton numbers, stating that the total number of fermions is preserved. Thus 
by considering a system which is so large that no fermion enters or leaves it, we can 
arrange that all physical observables satisfy f|3. 16j) . For a many-particle observable 
satisfying (|3.16p . the expectation value is again expressed by a trace, 

(0)=Tr? f (P f O) . (3.17) 

4. Microscopic Mixing of Decoherent Space- Time Regions 

In this chapter we develop a method for describing entangled states by a projector 
in the one-particle Hilbert space. In Section 14.11 we give a preliminary construction, 
which clarifies the difficulty in describing entangled states. In Section ^. 2l we overcome 
this difficulty on a rather formal and axiomatic level. The microscopic justification 
of the resulting formalism will be given in Sections I4.3H4.5I in the framework of the 
fermionic projector. Appendix IA.1I explains an attempt to describe entanglement by 
restricting attention to a subsystem. Seeing why this attempt fails might serve to the 
reader as a motivation for the concept of microscopic mixing. 

4.1. Microscopic Mixing of the Wave Functions. Our idea for realizing entangled 
states is to give P a non-trivial microscopic structure, with the hope that "averaging" 
this microstructure over macroscopic regions of space-time will give rise to an effective 
kernel P(x, y) of a more general form which allows for the description of entanglement. 
As the relativistic generalization will not be quite straightforward, we begin for clarity 
in the non-relativistic setting. Thus we consider the situation where space is subdivided 
into sets Mi , . . . , Ml , 

M = MiU---UM L and M a n M b = if a b, (4.1) 

which are fine-grained in the sense that every macroscopic region of space-time in- 
tersects several of the sets M a . The sets M a may be localized, but they can also be 
extended over a macroscopic region of space-time, for example by forming "layers" 
or "filaments" connecting the two observers in the spin correlation experiment of Ex- 
ample [34] (see Figure H]). The macroscopic physical objects are then introduced by 
homogenizing over the sets M a . We refer to this technique as the method of microscopic 
mixing. 

The partition (14. ID allows us to decompose "K into an orthogonal direct sum of the 
Hilbert spaces !K a of square integrable wave functions in M a , 

L 

M = 0M a . (4.2) 

a=l 
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Figure 1. Example of microscopic mixing in a spin correlation experiment. 



Thus every wave function in the image of P, (|3.8p . has the unique decomposition 



(a) 



(4.3) 



a=l 



The simplest attempt is to realize a macroscopic local one-particle observable O (like 
a position or spin operator) by an operator on "K being invariant on "K a , 



O : "K a — > "K a ■ 



(4.4) 



Then the corresponding one-particle expectation values split into a sum over the sub- 
systems, 

L L 

(iP\OiP) = V(V> (o) |Cty (o) ) =J2 i>{x)0^{x) dx , 

a=l a=l J Ma 

and this can be understood as an "averaging process" over the subregions M a . 

Following the constructions in Section 13.31 every one-particle operator induces a 
corresponding operator on the Fock space 3", and products of such operators yield 
corresponding many-particle observables. Taking expectation values of such operators 
in the Fock space again involves an "averaging process" over the subregions M a . More 
specifically, describing the many-particle system by a projector P in "K, the expectation 
value of the two-particle observables corresponding to Alice and Bob is given by (see 
Lemma I3.3() 

i=l o=l 
i,j=l a,b=l 



.,(&)> 



(4.6) 



In (|4.6p . an "averaging process" takes place at each of the observers. However, it 
is important to note that there is no averaging over correlations between the two 
observers, as would be the case for an expression like 



EE<#i°^! a) )(^i^ 

i,j=l a=l 



(4.7) 
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This shortcoming is the basic reason why the above method does not make it possible 
to realize general entangled states, as is worked out in Appendix IA.21 Thus our task 
is to extend our framework so as to obtain averages over correlations (|4.7j) . 

4.2. A Formalism for the Description of Entanglement. The previous construc- 
tion did not take into account that the measurement process is a result of an interaction 
of the system with the measurement device. Assuming that the subsystems have an 
independent dynamics (an assumption which will be justified in Section 14.41 below), 
also the measurement process should take place independently in the subsystems. Fol- 
lowing this idea makes it possible to describe entanglement, as we now explain. 

For the one-particle observables, the assumption of an independent dynamics of 
the subsystems was already taken into account in (|4.4|) by the assumption that O 
should be invariant on the subspaces "K a . However, for a many-particle observable, 
it was too naive to simply take the product of the one-particle operators (see (|4.5p 
and (14, 6h ) . Thinking of the situation in Figure [TJ Alice is built up of fermionic wave 
functions. Thus considering her as part of the physical system, we should replace the 
corresponding measurement operator Oa by separate operators for each of the 
subsystems. Proceeding similarly for Bob, in the subsystem M a measurements are 
to be carried out with the operators OjP and Og . For a correlation measurement 
in M a , we should extend the one-particle observables to operators 0^ a) defined 
on the Fock space 9^ of the subsystem given by 

3*0 = <A„(IK£)> c 3" n (4.8) 

(as explained after ()3.16|) . we again restrict attention to observables which preserve 
the number of particles) and consider the corresponding two-particle observable 

More generally, we make the following assumption: 

(A) The observables correspond to operators O which are invariant on 

O-.J^^J^, a = l,...,L. 

Similar as explained in Section [3?2l for the Fock space 3t , we can get a simple connection 
between a projector P in "K and the Fock spaces 9y • Namely, choosing again an 
orthonormal basis t/j\, . . . ,ipf of P(!K) and decomposing each of the one-particle wave 
functions according to (|4.3p . we can construct Hartree-Fock states in the /-particle 
Fock spaces of the subsystems, 

: = # A-A^fe 3j a) . (4.9) 

Exactly as in (|3.1U|) . one sees that these vectors are unique up to a common phase, 

w ith ip g R independent of a . (4.10) 

The setting so far is not sufficient for determining the expectation value of a mea- 
surement, because for computing an expectation value we need to take an "average" 
over the subsystems. This process can be described conveniently by the so-called mea- 
surement scalar product (.|.), which we now introduce (for a microscopic derivation of 
the measurement scalar product and the measurement process see Section [4.51 below) . 
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It is defined on the one-particle Hilbert space (3t, (•!•)) as a positive semi-definite inner 
product 

(.|.) : JfxM^C, (4.11) 

with respect to which the direct sum decomposition (|4,2p need not be orthogonal. The 
fact that this inner product is only semi-definite models the fact that the measure- 
ment process may involve a homogenization process on the microscopic scale, so that 
fluctuations of the wave functions on the small scale might not enter the measurement 
process. The measurement scalar product induces on the Fock spaces the bilinear form 

(.\.)^ : 3*0 X 3*0 -> C : (4 a) A • • • A ^\ # A • • • A 



^E(- 1 ) slgn(CT) (<l)i#)---(^ti^ ) : 



{O) = /t J . (4.12) 



n\ 

We now specify how expectation values are to be computed and state the assumptions 
which ensure that these expectation values are real. 

(B) The expectation value of the measurement of the observable O is given by 

Y, L b-i(* (a) 1 c ^W)( a ' b ) 

Ea, b =l(^ a) \^ b) Y a ' b) 

(C) The observables are symmetric possibly up to a microscopic error, meaning 
that 

(^(o) | 0& b) )? = (<D*W | *( fc ))gr + 0(e) , 
where e is the length scale of microscopic mixing. 
Finally, we need to specify how a state changes in a measurement process. In order 
to ensure that a repeated measurement of the same observable yields identical results, 
one usually asserts that after the measurement, the state should be an eigenstate of the 
observable. In our setting, the situation is a bit more involved because the measure- 
ment process may change the number of subsystems, and only the wave function after 
the homogenization should be an eigenstate of the observable. This is made precise 
by the following construction. We take the direct sum of the vector spaces 3"!^ , 

L 
a=l 

and on these spaces we introduce the inner product 

L 

(.|.)g:SxS-»C : ((* (a) )a=l,...,i,(* (6) )a=i,...,L)^ £(*(<0|tf (*))(«.*). 

a,b=l 

This inner product is positive semi-definite, but it need not be definite. Dividing out 
the null space and taking the completion, 

3f ■= STSo where % = {u G S with (u|n)g = 0} , (4.13) 

we obtain a Hilbert space, which we can regard as the effective n-particle Fock space 
obtained by homogenization over the subsystems. We denote the natural projection 
operator by ir n , 

L 

vrn : 03l a) ^?f • (4-14) 
a=l 
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Using linearity together with Assumption (C) above, every observable O induces an 
operator 

O cS : 3f ->3f , (4.15) 

uniquely defined possibly up to a microscopic error. To the projector P in JC we 
associate a corresponding state 

q> cS = 7 r / (V 1 ),...,¥ i )) G 3"f (4.16) 

(with the VPW as in (|4.9p ). According to (|4.10p . this state is well-defined up to an 
irrelevant phase. 

(D) After a measurement of the observable O, the one-particle projector P takes 
such a form that the corresponding state ^ G 3~j ff defined by (|4.16p is an 
eigenstate of the operator eff , (|4.15p . 

Similar as in the Copenhagen interpretation or the formulation of the measurement 
process by von Neumann [20j, Section V.l], the above Assumptions (A)-(D) are merely 
working rules to determine the results of measurements. For a conceptually convincing 
treatment, these assumptions should be derived from the physical equations. 

We now verify that the above setting indeed makes it possible to realize the EPR 
singlet state. 

Example 4.1. (The spatially separated singlet state) We choose a microscopic length 
scale e > and subdivide position space M = R 3 into two subregions M\ and M2 
which form layers of width e, 

Mi ={i£l 3 with [xx/e] G 2Z} , M 2 = {fef 3 with [x x /e] G 2Z + 1} 

(where [x] = min{n G Z | n > x} is the Gauss bracket). We introduce the wave 
functions 

V>i(x) = ip\(x) XMi(x) + ip A (x) XM 2 (x) 

I "t" 

ip2{x) = ijr B {x) XM x {x) -ip B (x) xm 2 [x) , 

where ip^ B are smooth one-particle wave functions supported near Alice or Bob (and 
where x is the characteristic function defined by Xn{%) = 1 if x G N and Xn{%) = 
otherwise). Defining P as the projector on the subspace spanned by tp\ and ip2, the 
corresponding two-particle wave functions of the subsystems are 



(1) = c V^A A ^b) XMixMi and * (2) = -c A 1/)%} XM 2 xM 2 



with c a normalization constant. 

In order to realize a suitable mixing of the subregions in the measurement process, 
we introduce the measurement scalar product by 

(ip\(j>) = f Tp(x)(f)(x) dx + — I (ip{x + ee\)(j)(x) + ij){x)<j>{x -\- £e\)\ dx . 
Jm 2 J M \ J 

The spin operators are symmetric with respect to this inner product, whereas the 
position operators are symmetric up to an error of order e, 



xip\<j)) — (ip\x(f>) = - xi (ip(x + ee\)(j)(x) — ip(x)<f>(x + ee.\) ) dx . 

2 J M \ ) 
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Thus the general observables introduced according to (A) indeed satisfy the condi- 
tion (C). The expectation values of the spin operators can now be calculated by ap- 
plying the rule (B). More precisely, the inner products in (|4.12p are computed by 

(^(l)|0$(l))(i.i) 

(^(2)|0$(2))(2,2) 
($(1)|0$(2))(1,2) 

(^(2)|0#(1))(2,1) 

where on the right the scalar product on J as defined by (|3.ip appears, and the 
subscript + denotes that both spatial arguments of the corresponding two-particle 
wave function have been shifted by ee±. Note that all these inner products involve 
integrals over Mi x Mi or M 2 x M 2 . Since the wave functions ip^/ B are all smooth, 
we can extend the two-particle wave functions of the subsystems to smooth functions 
on M x M, 

*2 = c(v>Wi) and *2 = c (ri A 4)- ( 4 - 17 ) 

Shifting the arguments changes these smooth wave functions only by a term of order e. 
Also, extending the integration range in the above integrals from M\ x M\ or M 2 x M 2 
to M x M changes the values of the integrals only by a factor of four, again up to 
contributions of order e. We thus obtain 

L 

£ («) I O = 4 + I O + *®)) + • 

a,&=l 

We conclude that in the limit e \ 0, the expectation values as defined by (B) indeed 
coincide with the expectation values of the spin singlet state. 

Moreover, it is straightforward to verify that the space 3"| ff as defined by (|4.13f) 
can be identified with the ordinary Fock space 3~2 ; an d that under this identification, 
the state ^ efT as given by (|4.16p goes over to the state + w ^h the wave 
functions as in ()4TT7D . 

4.3. Microscopic Mixing in the Framework of the Fermionic Projector. We 

now generalize and adapt the method of microscopic mixing to the framework of 
the fermionic projector. For the extension to the relativistic setting, we decompose 
Minkowski space M into two disjoint subregions 

M = M 1 UM 2 with M 1 DM 2 = 0, (4.18) 

which may be fine-grained as depicted in Figure [2] (for simplicity, we only consider two 
subregions; the generalization to a finite number of subsystems is straightforward). 
If one prefers, one can in addition replace the space-time continuum by a discrete 
set of points, but the distinction between continuum and discrete space-time will not 
be of relevance for the following considerations. We again consider a family of wave 
functions ip\, . . . ,4>f, where in view of the fact that we also count the states of the 
Dirac sea, the number / of particles is very large (see Section I2.3|) . As in ()2.5p , the 



= (*( 2 )|£>$(2)) 

= i((*2 ) |0*( a )) + (*W|0*2 ) )) 
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Figure 2. Example of a microscopic mixing of two space-time regions. 



fermionic projector takes the form 

/ 

P(x,y) = -y2\ipj(x)y^.ipj(y)\ with x,yeM 1 UM 2 . (4.19) 
i=i 

Splitting up the wave functions similar to (|4.3p by 

^• = ^ 1) +Vf ) with ^f=^ jXMa (a = 1,2), (4.20) 

to each subsystem we can associate similar to fj4.9j) the many-particle wave function 

M = ^ A . . . A . (4.21) 

At this point, we need to discuss the significance of the inner products (|2.8p and 
(I2.10p . In discrete space-time, the underlying inner product involves the sum over all 
space-time points (|2.3p . Likewise, in the continuum limit, the sum is to be replaced by 
a space-time integral (|2.8p . The notions of symmetry and idempotence of the fermionic 
projector refer to this indefinite inner product (for details see O §2.6]). However, as 
the fermionic projector is built up of negative definite states, we may always restrict 
attention to a negative definite subspace, on which (.|.) := — <.|.> is a scalar product. 
After completion, we again get a separable Hilbert space (5f, (.|.)), so that we are back 
in the setting of Chapter [3j It is important to keep in mind that the inner product (|2.8p 
is not suitable for computing the expectation value of a measurement, which usually 
takes place at a fixed time. Thus in the measurement process, it is natural to work 
with a different scalar product, which can be regarded as the relativistic analog of the 
measurement scalar product (|4.1ip . In view of the fact that the integrand of (12. 10p 
has the interpretation as the probability density, the scalar product (I2.10p seems the 
correct choice. 

Let us consider how to implement the scalar product (12. 10f) as the measurement 
scalar product in the presence of a microscopic mixing of two subsystems. First, 
taking into account that realistic measurements take place in a finite time interval, it 
seems a good idea to replace the spatial integral in (|4.18p by an integral over a strip of 
width At (as shown in Figure[2|). Moreover, the measurement should involve some kind 
of homogenization process on the microscopic scale. This is clear empirically because 
the expectation value must be a computable quantity which involves taking averages 
over the subsystems. The homogenization can also be understood microscopically from 
the fact that the measurement devices are themselves formed of quantum mechanical 
wave functions which are spread out in space-time. A typical example for an idealized 
measurement device is the operator \rj> <rj\, where n is a wave function which is 
supported inside a strip of width At as shown in Figure [2j The simplest way to take 
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into account the effect of such a measurement device would be to consider instead 
of (|2.10p the measurement scalar product 

i r rt+At rt+At 

^ = 7X7^2 M dt ± dt 2 ^(t 1 ,x)\-f cf ) (t2,x)y , (4.22) 

{Aty j R3 j t j t 

where we take the average over a small time interval before computing the spatial inner 
product. A more realistic measurement device is of course much more complicated, 
but fortunately the details are of no relevance. All that matters is that the measure- 
ment scalar product involves a homogenization process, with the effect that the direct 
summands in (14. 2D are in general not orthogonal with respect to (.|.) and that the null 
space So m (14.13P becomes non-trivial. 



4.4. Justification of an Independent Dynamics of the Subsystems. Following 
the concept of discrete space-time of Section 12. 1\ space-time is not smooth on the 
microscopic scale, but should have a non-trivial microstructure. In order to explain 
the possible consequences of such a microstructure in a simple setting, we consider what 
happens if we choose the wave functions ip^ in the two subregions differently. More 
specifically, we transform the wave functions in the second subsystem by a unitary 
matrix of determinant one, 

/ 

Vf } -> , -> $ 2) := U Jk V>f with U G SU(/) . (4.23) 

k=l 

This transformation has the advantage that it has no effect on the many-particle wave 
functions (I4.2ip . because (exactly as in (I3.10j) ) 

$(2) §(2) = ^(2) A . . . A ^(2) = ^ jj ^(2) A . . . A ^2) = ^(2) _ 

Furthermore, the fermionic projector does not change if its two arguments are in the 
same subsystem, because in the case x, y G M 2 , the unitarity of U yields that 

/ / 

j=l j,k=l 

f 

= _J2\4,f(x)y^f(y)\=P(x,y). 
j 

However, if the two arguments of the fermionic projector are in different space-time 
regions, the operator U does not drop out. For example, if x G M 2 and y G Mi, 

/ / 
P(x,y) ^-Y^\^f\ x )^f\y)\ = - U jk \il>f\x)^£\v)\ ■ (4-25) 

j=l j,k=l 

In the special case when U is a diagonal matrix whose entries are phase factors, 

/ 

U = diag(e^ 1 , . . . , e iip f ) with ^ ^ = mod 2ir , 
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the summations in (|4.25p reduce to one sum involving the phase factors, 

/ 

P(x,y) -> -J^e^ \^ 2) (x)y^f\y)\ . 

If the angles ipj are chosen stochastically, the phases of the summands are random. 
As a consequence, there will be cancellations in the sum, and keeping in mind that 
the number of summands is very large, we conclude that P(x, y) will be very small. 
More generally, we find that if U is a random matrix, P(x, y) becomes small if x and y 
lie in different subsystems (this argument will be quantified in Section 14.61 below by 
integrating over the space of unitary matrices). 

From the physical point of view, the above consideration can be understood using the 
notion of decoherence. If the one-particle wave functions ipj and ipj are coherent 
or "in phase", then the fermionic projector P{x,y) has the usual form, no matter 
whether x and y are in the same subsystem or not. If however the wave functions in 
the subregions are decoherent or "out of phase", then the fermionic projector P{x,y) 
will be very small if x and y are in different subregion. We refer to this effect as 
the decoherence between space-time regions. It should be carefully distinguished from 
the decoherence of the many-particle wave function (see for example [Hj). Namely, 
as we saw in (|4.24[) . in our case the many-particle wave functions remain unchanged. 
Thus they remain coherent, no decoherence in the sense of [13] appears. But the one- 
particle wave functions become decoherent (|4.23p . having an influence on the fermionic 
projector (|4.19p . 

We next consider the influence of the decoherence between space-time regions on 
the dynamics of our system. We begin by discussing the extreme case where P(x, y) 
vanishes identically for x and y in different subsystems. Then the action (|2.7p splits 
into the sum of the actions of the two subsystems, so that the interaction takes place 
independently in the two subsystems. In other words, the subsystems decouple. By 
restricting two different systems in Minkowski space to Mi respectively M2, one can 
apply the methods of Section [231 to both subsystems separately. Then each subsystem 
is described by an independent continuum limit in terms of a Dirac equation (|2.11|) 
coupled to a classical field (|2.13p . This explains the assumption of an independent 
dynamics of the subsystems made in Section 14.21 We point out that, following the 
concept that the wave functions generate the causal and geometric structure of space- 
time (see the last paragraph in Section 12. 3p , the decoupling of the subsystems even 
implies that between the subsystems, the usual causal and topological structure of 
Minkowski space ceases to exist. 

If we merely assume that P(x, y) is small for x and y in different subsystems, our 
action principle (|2.7p does describe a coupling of the two subsystems, which however 
should be weak. Keeping in mind that the causal structure of Minkowski space is 
related to the singularities of distributions like (12. 9p on the light cone, and that this 
singular structure will be destroyed by decoherence, we know that the coupling of 
the two subsystems cannot be described by causal equations formulated in Minkowski 
space. Although we have a precise mathematical framework (12. 7p . describing the 
coupling of the subsystems quantitatively seems very difficult and goes beyond the 
scope of this paper. But from the mathematical structure of our action it is already 
clear that we do not get contributions from the boundaries of the two subregions. 
Therefore, instead of considering the "layers" in Figure [IJ it seems more appropriate 
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to draw each subsystem as many disconnected "bubbles" in space-time as shown in 
Figure [2j In view of the continuum limit, each system has an underlying smooth 
structure inherited from a corresponding system in Minkowski space, as is indicated 
in Figure [2] by the two "smooth space-time sheets" Afj^ nt . But of course, this picture 
should be considered only as a vain attempt to illustrate an unknown and probably 
very complicated microscopic structure of space-time (a maybe more realistic picture 
will be outlined in Section |5\3|) . 

4.5. Justification of the Superposition of Fock States. We now want to verify 
that expectation values computed with respect to the measurement scalar product 
indeed involve superpositions of Fock states. For simplicity, we again consider the sit- 
uation for two subsystems Mi and and assume that the fermions in each subsystem 
are described by an n-particle Hartree-Fock state 

whereas the remaining f — n particles describe the Dirac sea. Thus we choose the 
one-particle wave functions before microscopic mixing as ip[ , . . . , ipn , V'i+l' • • • > V*/ > 
where the first n wave functions describe the particles, whereas the other wave func- 
tions form the sea. We again introduce a decoherence between the subsystems by a 
unitary transformation of all states in the second subsystem (|4,23p . 

Expectation values (ipi\ipj) of the one-particle wave functions with respect to the 
measurement scalar product involve a homogenization process, with the result that 
wave functions should be identified which differ only by microscopic fluctuations. More 
specifically, we should not distinguish between the sea states of the two subsystems. 
Thus introducing on "K the equivalence relation 

ijj ^ $ - $ | ip - tP) = , 

we assume that 

^(i)^2) for allj = n+l,...,/. 
For ease in notation, we make this identification clear simply by omitting the corre- 
sponding superscripts W or ^ . 

Under this identification, the many-particle wave function of the whole system be- 
comes 

^ = tpi A • • • A tpf 

= + ^i 2) ) A • • • A (V« + V^ 2) ) A (Vn+1 + Vwi) A • • • A (ipf + ipf) . (4.26) 

Multiplying out, we obtain many contributions. One of them corresponds to the many- 
particle wave function of the first subsystem 

4 1] A • • • A V4 1} A Vn+i A • • • A V/ , (4.27) 
and another is the many-particle wave function of the second subsystem 

$ 2) A • • • A tp^ A 4> n+1 A • • • A i) f = det U ip^ A • • • A ip^> A ip n+ i A • • • A ip f 

= 4 2) A • • • A A 4> n+1 A • • • A V/ • (4.28) 

All the other contributions involve matrix elements of the unitary operator U. Similar 
as explained after (14. 25ft , all these contributions become small if U is a random matrix. 

We conclude that the measurement process involves the sum of the many-particle 
wave functions (|4.27p and (|4.28p of the two subsystems. This justifies Assumption (B) 
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in Section \4. 21 Moreover, this consideration explains why for measurements one should 
work in the Fock space 3f defined by (ETT31) . 

4.6. Describing Fock Superpositions with Random Matrices. In order to make 
the consideration of the previous section more precise, we now reformulate its mathe- 
matical core in terms of random matrices. Note that SU(/) is a compact Lie group, on 
which we consider a probability measure d\i. Then taking the average over a random 
matrix U G SU(/) corresponds to integrating U with respect to d[i. The simplest 
choice for d\i is the normalized Haar measure d^Haar (see for example [2J Section 1.5]), 
but other choices are possible, as will be discussed below. 

We first observe that certain products of matrix elements of U vanish on average. 
For simplicity, we consider the Haar measure. 

Lemma 4.2. Suppose that for any p in the range 1 < p < / — 1, we choose in- 
dices ii,..., i p and j±, ■ ■ ■ ,j p with i\ < ■ ■ ■ < i p . Then 

SU(/) 

Proof. We let k be an index which is not contained in . . ,i p } and let V be the 
diagonal matrix which has entries one, except that = V^k = — 1- Then V 6 SU(/), 
and thus a variable transformation shows that the above integral is invariant under 
the replacement U — > VU. But this transformation flips the sign of the integrand. □ 

Applying this lemma to the expression (|4.25p . we see that the the fermionic projec- 
tor P(x, y) indeed vanishes for x and y in different subregions, if the mean value 
over SU(/) is taken. The lemma also applies to the contributions obtained by mul- 
tiplying out (]4.26p . It shows that all contributions vanish on average, except for the 
many-particle wave functions (I4.27P and (I4.28P of the two subsystems. 

Since the expectation value of a measurement involves the absolute square of the 
wave functions, we also need to integrate the absolute square of the many-particle wave 
function (j4.26j) over SU(/). We begin with a simple integral involving the absolute 
square of one matrix element of U. 

Lemma 4.3. For any j, k £ {1, . . . , /}, 

,^fc| 2 dflHaar = 7 • (4.29) 
SU(/) / 

Proof. By multiplying with suitable unitary operators from the left or the right, we 
can arbitrarily change the values of the indices j and k, without changing the integral. 
Thus 

/ 



/ \Ujk\ 2 ^Haar = 77 / Wjk\ 2 C^Haar 

■/sue/) J J sv(f) ~k~ti 

= 72 I Tr(U*U) d^Haar = 7 
/ JSV(f) J 



'SU(/) 

concluding the proof. □ 

Applying this result to (|4.25p , we see that decoherence typically scales the kernel of the 
fermionic projector down by a factor f~2 . This quantifies that P(x,y) really becomes 
small if x and y lie in different subregions. 
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Lemma 14.31 could be generalized to integrals involving the absolute squares of n 
matrix elements, giving the result 

f , ,2 , ,2 1 

/ Fiifc! I • • • \Uj n k n \ d/i H aar ~ -JZ if fl < / . 

This shows that every fixed summand obtained by multiplying out (|4.26p except 
for (|4.27p and (|4.28p vanishes in the limit / — > oo. Unfortunately, this is not quite 
good enough, because the number of summands becomes large if / increases. Thus in 
order to estimate the whole sum of terms, we need to use a different method, which 
we now explain. 

To describe the combinatorics of the wave functions tp^ and ifi\ 2 \ we consider a 
subset jW C {1, . . . , n} and take its complement 1^ = {1, . . . ,n} \ /W. As the case 
of no particles is trivial, we may assume that n > 1. Suppose that we are interested 
in the contribution to (|4.26p of the form 

* - < u ) ( A i6 *) A ( A je/(2 , 4 ] ) A A ••• A ^ 

with a complex prefactor c(U). Multiplying out only the first n factors in (|4.26|) and 
using the definition of fo 2 ' in ()4.23p . we find that 

c(U) = sign(/«) det(X^ + UX&), 

where are the diagonal matrices 

(x^)* =^ X/(a)u{n+w} « 

(here x 1S the characteristic function and sign(J) defined by 

sign(J) = (_l)«i+"-H»+ a£ T li 

if we consider / as the ordered set I = (ft, . . . , i g ) with 1 < i% < %i < ■ ■ ■ < i g < /). 
Thus to take the average of ^ and \^\ 2 , we need to compute the integrals 

2 



/ det(X (1) +UX {2) ) d/i and / det(X (1) + C/X (2) ) 



dp, , 



'su(/) Jsv(f) 

respectively. We get agreement with the Fock space formalism if and only if the 
following identities hold: 



lim / det(X« + UX®) dfi = 5 lW + 5 jW (4.30) 

f-*°°Jsv(j) 

lim f det(XW + UXW) det(X (1 ^ + UX®) dp 



c $imjm [ s im,0 + <W ) • (4.31 



Here the parameter c > is an overall constant (to be chosen independent of /), which 
can be absorbed into the definition of the scalar product on the Fock space 9\^ ff . 

In the above integrals we worked with a general probability measure d[i. Let us 
discuss how it is to be chosen. Ultimately, this measure should be determined by 
analyzing the statistics of the decoherent subsystems which form according to the 
action principle (|2.7p . Since such an analysis is not yet available, we must rely on 
heuristic considerations. One guide line is that d\x should respect all symmetries of 
the underlying framework. The simplest measure which meets this requirement is the 
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Haar measure d^Haar- Writing the above determinant as de^lA^ 1 ) + UX^), one sees 
that the identity matrix is distinguished, and thus we can choose dfi more generally 
as any measure which is formed of U and 1. For example, one could choose dfi to be 
a constant times the measure 

| Tr(l + U)\ a d/XHaar, |Tr((l + U) a ) | d/i H aar Or | det(l + U)\ a d/i H aar 

with a£l, but many other choices are possible (see for example |16|). The question 
for which choices of dfi the relations (|4.30[) and (|4.3ip hold is an open problem which 
will be considered elsewhere. 

We finally remark that the above considerations immediately generalize to more 
than two subsystems, if one replaces the term det(X^ + UX^) by 

det(xW + U 2 xW+--- + U L xM) 

and integrates over all random matrices U2, . . . ,Ul G SU(/). 



5. Second Quantization of the Bosonic Field 

In Section 14.41 we considered two decoherent subsystems M\ and M2 and saw that 
by analyzing each subsystem in the continuum limit, we could describe the dynamics 
by the Dirac equation coupled to a classical field. Taking a finite number of such 
decoherent subsystems, the whole dynamics is described by several classical fields, one 
for each subsystem. In this chapter we will show that the resulting framework indeed 
allows for the description of second quantized bosonic fields. 

For simplicity, we merely consider an electromagnetic field (the generalization to 
other bosonic fields is straightforward). We subdivide Minkowski space into L disjoint 
regions M\, . . . ,Ml, which are again assumed to be fine-grained. Similar to ()4.19j) 
and (|4.20p . the fermionic projector can be written as 

/ L 
P(x,y) = -J2\^>^j(v)\ with ^ = 5>j a \ ^ a) -=^XM a , (5.1) 

j=l a=l 

where / is a large number which tends to infinity if the ultraviolet regularization is 
removed. As in Section \AA\ we arrange by unitary transformations of the form (I4.23D 
that the subsystems are decoherent. Considering each subsystem in the continuum 
limit, we obtain similar to (12. lip and (I2.13P the Dirac-Maxwell equations 



(iyi(dj - ieA { f) - rnjP {a \x,y) = 



n f n a (5.2) 

k=l 1=1 

We note for clarity that according to (|4.20p . the wave functions and <j>f^ are 
obtained by restriction to a subregion M a C M of space-time. But it is reasonable 
to assume that these wave functions are macroscopic in the sense that they can be 
extended smoothly to the whole Minkowski space. Similarly, we assume that the 
potentials and the fermionic projectors are defined on a whole sheet M^ ont 
of Minkowski space (see Figure [2|) . 
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5.1. Describing a Second Quantized Free Bosonic Field. In order to get a 
simple connection to standard textbooks like |23[ 1 15 j . we begin with a free electro- 
magnetic field (i.e. the situation where no fermionic particles or anti-particles are 
present). Furthermore, to avoid the technical issues involved in taking an infinite vol- 
ume limit, we restrict attention to the situation in finite spatial volume by considering 
a three-dimensional box of length I with periodic boundary conditions. Working in 
the Coulomb gauge divvl = 0, the Maxwell equations reduce to the ordinary wave 
equation for each component of the vector potential, 

DA(t, x) = 0, 

whereas the electric potential A can be arranged to vanish identically. Decomposing A 
into Fourier modes of momentum k £ (2-7rZ/£) 3 , 

A(t, x) = Y, k) e lfe + a(t, k) e" ife ) , 

k 

the Maxwell equations are solved by 

a(t, k) = a(k) e~ iujt with co := \k\ , 

whereas the Coulomb gauge gives rise to the transversality condition k ■ a{k) = 
(see [15, Chapter I, §2]). The two linearly independent solutions of this transversality 
condition correspond to the two polarizations of the electromagnetic wave; we denote 
them by an index (3 = 1,2. Introducing the canonical field variables 

= 4"- ( a ^) + a /?(>)) > P(*) = ^0$) = (ap(k) ~ OflW) , 

the energy H of the classical electromagnetic field becomes (for details see |15l Chap- 
ter I, §2]) 

H = E £ \{v P (k? + ^(k?). (5.3) 

fce(27rZ/£)3 0=1)2 

Here each summand is the Hamiltonian of a harmonic oscillator. Thus we have rewrit- 
ten the classical electromagnetic field as an infinite collection of classical harmonic 
oscillators. 

The second quantization of the electromagnetic field corresponds to quantizing each 
harmonic oscillator as in standard quantum mechanics (see for example |23t Part I, 
Section 1.2]). We proceed by discussing the connection between the classical and 
the quantum dynamics in detail, for simplicity for a single harmonic oscillator of fre- 
quency uj. Thus our starting point is the classical Hamiltonian 

h(p,q) = ^{p 2 +u;V). (5.4) 

Here q and p are the canonical variables, which together form the classical phase 
space V = {(p,q) with p,q E R}. The classical dynamics is described by Hamilton's 
equations 

d dh d dh 2 

dt^ dp ^ ' dt^ dq u 1 ■ 
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A solution (p(t),q(t)) describes a classical trajectory. Solving Hamilton's equations, 
the classical dynamics describes a rotation in phase space, 

( P(t) \ _ (cos cat — sinwA / p(0) 
\ojq(t) J ysinut cos ujt J \uq(0) 

In order to get a setting similar to that in quantum theory, we next consider on 
phase space complex- valued functions ip(p, q), referred to as "classical wave functions." 
Introducing the scalar product 



(5.5) 



(i>\4>) class = // i>(p,q)<t>(p,q) dpdq , (5.6) 



the classical wave functions form a Hilbert space (% c iass> (■lOciass)' The phase flow f|5.5p 
induces a flow on H, which is described most conveniently by the time evolution oper- 
ator Udass defined by 

EWt)^)(p(t),?(*)) = (UciassiQ) ^)(p(0),g(0)) • (5.7) 

It is a unitary operator on % c i ass . Before going on, we remark that in classical physics 
one usually works instead of complex functions with positive functions or densities 
on phase space. Working with complex-valued functions and the scalar product (15.61) 
seems unusual but will be very useful for the following considerations. In a some- 
what different context, the Hilbert space (% c iass> class) is also used in geometric 
quantization for the the so-called prequantization (see [22], Section 5.2]). 

For the quantization of the oscillator, one replaces the canonical variables p and q 
by self-adjoint operators P and Q which act on a complex Hilbert space (H, {■{■)) 
and satisfy the canonical commutation relations [P, Q] = —i. The physical system is 
now characterized by a state £ T~L. The dynamics is described by the Schrodinger 
equation 

id t ty = with H = - (P 2 + u 2 Q 2 ) . (5.8) 

It is most common to represent T~L as the space of square integrable functions with the 
inner product 

(tf |$) = / dq , (5.9) 

Jr 

and to choose the operators Q and P as 

m)(q)=q^(q) and P ' = . 

dq 

Integrating the Schrodinger equation gives rise to the unitary time evolution operator 

U(t) = e~ itH : Oi -> Ji : *(0, q) ^ q) . (5.10) 

With the above formulation we expressed both the classical and the quantum dy- 
namics by a unitary time evolution operator acting on a Hilbert space (see (|5.7p . (j5.6[) 
and ()5.10p . (|5.9fl ) . But the time evolution operators have a completely different form. 
Furthermore, the Hilbert spaces are different, because the "classical wave functions" 
depend on both q and p. In the quantized theory, however, the Heisenberg uncertainty 
principle prevents P and Q from being simultaneously measurable, as is reflected 
mathematically by the fact that they correspond to non-commuting operators. Since 
in the classical theory, position and momentum can be chosen independently, there is 
much more freedom to choose the initial wave function ip(p, q) than in quantum the- 
ory, where choosing ^(q) automatically determines the corresponding wave function in 
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momentum space. This raises the question if for a given quantum wave function ^(q) 
we can choose a corresponding classical wave function ip(p, q) such that the classical 
dynamics of i/j as described by (|5.7p coincides with the time evolution of the quantum 
wave function (|5,10p . While the general answer to this question is no, it turns out that 
for the harmonic oscillator this correspondence can indeed be made: 

Lemma 5.1. (Correspondence between classical and quantum dynamics) 

Consider the classical harmonic oscillator (|5,4p with dynamics (|5.5p , (|5.7p and the cor- 
responding quantum harmonic oscillator with the dynamics described by the Schrddinger 
equation (|5,8p and (|5.10p . Then there is an isometric embedding i : % — > H c iass which 
maps the quantum evolution to a corresponding classical evolution, in the sense that 

U c ia SS {t) l = iU{t) e iwt/2 forallteR. (5.11) 

Moreover, there are differential operators H c i ass , P c i ass and Q class in H class such that 

H class L — ^ H , Pciass ^ — ^ P and Q class ^ — I Q • 

We point out that the factor e lujt / 2 in ()5.11|) corresponds to the zero point energy of 
the quantum harmonic oscillator. This factor modifies the wave functions only by a 
joint global phase, without an influence on any observations or expectation values. 

Proof of Lemma \5.1[ We choose an orthonormal eigenvector basis of the Hamil- 
tonian in ()5.8[) (see for example [181 Section 3.1]) 

HV n = (n + ~ J w* n , n = 0,l,.... 

Writing the Hamiltonian as H = ui{a)a + ^) with the annihilation and creation oper- 
ators 

' r wq+4-)> at = ^f^-#V (5.12) 



y/2u V dqj ^phfj V dq / 

the eigenvectors can be obtained by acting with the creation operators on the ground 
state, 

= c exp (~^f-) and = c n (a f ) n , (5.13) 

were the c n are positive normalization constants. From (|5.10p it follows immediately 
that 

U(t) e iu}t / 2 ^ n = e~ inujt * n . (5.14) 

In order to define the mapping l, it suffices to associate to every eigenfunction ^> n a 
corresponding classical wave functions ijj n G % c iass (then i is determined uniquely by 
linearity and continuity). First, in order to write the classical dynamics in a simpler 
form, we rescale the momentum variable by introducing the new phase space variables 

x = q and y = — . 

Setting z = x + iy, the time evolution operator (j5.7|) becomes 

(U c i ass m)(z) =i>{e wt z). (5.15) 
We now define the "classical annihilation and creation operators" on % c i ass by 

Oclass = ^ + ia v) > "class = \ (<4 + ia \ 
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where a x and at are given in analogy to (|5.12p by 



1 ( d\ f If d 

oox + — , al = , oox 



\/2uj V dx J ' x y/2uj V dx t 
whereas the subscript y refers similarly to the variable y. We introduce the wave 
functions ip n by 

2 ) and V'n = c n (aJ lass ) n ^o • (5.16) 

Let us verify that the resulting mapping t has the required properties. First, it is 
obvious from their definition (|5.16p that the functions tpn are orthonormal in % c i ass , 
and thus i is indeed an isometric embedding. Using a polar decomposition z = re tip , a 
short calculation shows that 



idfi a class 



class 



Applying this relation in (15.16P and using that is radially symmetric, we obtain 

id^t/jn = nip n and thus Vv^) = e~ inip (f> n (r) 
with radially symmetric functions <f> n . Thus the classical dynamics (|5.15p implies that 

U cl3SS (t) ^ n = e-^ nt 

Comparing with (|5.14p proves (|5.1ip . 

In order to construct the operators H c \ ass , P c i ass and Qdassi w e first note that both 
the classical and quantum annihilation and creation operators satisfy the canonical 
commutation relations 



a class> °class 



a\a 



1. 



1 and 

and in view of (|5. 13[) and (|5.16p they correspond to each other in the sense that 

a c iass t = i a and aj lass i = ia} . 

Thus expressing the operators in % in terms of a and a', we obtain the corresponding 
"classical" operators simply by adding subscripts. Thus we set 

#class = 00 ^ lass Oclass + - 

and 

Qclass = r— (a c l ass + a£ lass J , -Pciass = — ( a class _ a ^ ass 



/2a; V ciassy . — V 2 

This concludes the proof. □ 

We remark that the mapping i appears in the mathematical physics literature as the 
so-called Bargmann transform (see |2H Section 4.3]). But to our knowledge, it has not 
been used to get a connection between the classical and quantum dynamics. 

The above lemma shows that by choosing the "classical wave function" <j) E % c iass 
appropriately, we can arrange that the classical dynamics reproduces any quantum 
dynamics. In simpler terms, the quantum dynamics of the harmonic oscillator can 
be recovered as a special case of the classical dynamics. However, for making this 
correspondence, we had to take a somewhat unusual point of view and work on clas- 
sical phase space with complex- valued functions and the scalar product (|5.6p . To us, 
Lemma 15.11 is useful because it makes it possible to approximate a quantum state by 
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a finite number of classical trajectories, if to every classical trajectory we associate a 
complex number. This can be seen as follows: Suppose that we want to describe a 
quantum state ^ 6 "K. According to Lemma 15.11 this state has the same dynamics as 
the classical wave function t/j := ity G Ji c \ ass . For any L E N and an index o = 1, . . . , L, 
we now choose points (p^ , q^) in phase space together with complex coefficients 4>(a) 
which approximate if) in the sense that 

L 

<K°) KP - P (a) ) S (l - ? (a) ) ^ q) , (5-17) 

a=l 

with convergence in the distributional sense. For these discrete configurations, we can 
make sense of the scalar product (|5.6p by setting 

a,b 

Thus by choosing L sufficiently large, we can approximate the quantum dynamics of \& 
by a complex valued function (j) defined on a finite number of classical trajectories. 

Before going on, we point out that the scalar product (15.61) is invariant under the 
local phase transformations 

feg)^e^^( M ), (5.18) 

where cp is a real-valued function on phase space. Thus the phase of the functions 
in ^Kciass has no physical relevance. What counts is only the relative phase when 
taking superpositions of two wave functions ip,4> G % c iass- Similarly, in the discrete 
approximations in (I5.17p . the phase of the function <j>(a) can be changed by 

0(a) -> e tv{a) (t){a) under the constraint (p (a) ,q (a) ) = {p {b) ,q (b) ) =>- (p(a) = <p(b) . 

In the remainder of this section, we extend the above considerations on the harmonic 
oscillator to the Hamiltonian of the electromagnetic field (|5.3p . By taking tensor prod- 
ucts, the result of Lemma 15 . 1 1 immediately extends to a collection of harmonic oscilla- 
tors as in (|5.3|) . It then states that by considering suitable complex- valued functions 
on the set of all classical field configurations, one can reproduce the full dynamics of 
the free second-quantized field. Using an approximation argument similar to (|5.17|) . it 
suffices to consider a finite number of classical field configurations. Thus the remaining 
task is to associate to every classical field configuration a complex number (fi(a). Let 
us return to the setting of decoherent subsystems in the continuum limit (|5.2[) . Then 
every subsystem involves a classical electromagnetic potential A^ a \ In the considered 
case without fermions, the field equations reduce to the free Maxwell equations, i.e. 
again in the Coulomb gauge 

al ( - a \t,x) = (5.19) 

and (A^)°(t, x) = 0. Moreover, we have the Dirac equation for the fermionic projec- 
tor, which according to (I2.12D consists only of the sea states, 

(i~f j (dj - ieA { 3 a) ) - mj P sca (x, y) = if x e M a . (5.20) 

Having an ultraviolet regularization in mind, the number / of sea states is finite, so 
that P sea can be written in the form (15. ip . For a given solution A of the free Maxwell 
equations, we now introduce the following reference system. The causal perturbation 
expansion distinguishes a subspace of the solution space of the Dirac equation as being 
formed by the sea states. Selecting the /-dimensional subspace of the sea states which 
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is compatible with our regularization and choosing an orthonormal basis ipi , . . . , ipf 
of this subspace, we can introduce the many-particle wave function tp of our reference 
system by 

* = $1 A • • • A-0/ . 

As in (|3.1U|) . the freedom in choosing the orthonormal basis implies that is deter- 
mined only up to a phase. Now suppose that A coincides with the electromagnetic 
potential A a in one of our subsystems. Then the wave functions tp^ , . . . , ipi obtained 
by restricting the wave functions of the fermionic projector of the whole system to the 
subsystem M a span the same subspace of the Dirac solutions as the vectors ipi, . . . , ipf 
(probably after suitably modifying the solutions on the microscopic scale or modifying 
the regularization; a technical issue which for simplicity we will ignore here). Hence the 
corresponding many-particle wave function (|4.2ip coincides with W up to a complex 
number, 

if( a ) = (f>( a ) ^ with (f>(a) G C . (5.21) 

In this way, we have associated to the field configuration A a a complex number (p{a). 

Let us consider the phase freedom. As noted above, the phase of the wave func- 
tion depends on the choice of the basis ipi , . . . , ipf . Similarly, by transforming the 
orthonormal basis ipi, . . . , ipf of the image of P, we can also change the phase of 
arbitrarily. Thus (|5.2ip is well-defined only up to a phase. Now suppose that A also 
coincides with the electromagnetic potential A b of another subsystem. Then writing 
the many-particle wave function of the new subsystem as = cp(b) transforming 
the bases tp%, . . . , ipf or ip%,...,ipf changes the phase of both <p{a) and <p{b) in the 
same way. Thus the relative phase of (p(a) and (p{b) is well-defined. In other words, 
the complex- valued function <p is uniquely defined up to the transformations 

4>{a) -)■ e i¥,(a) (p(a) under the constraint A a = A b => ip(a) = ip{b) . (5.22) 

These transformations can be regarded as local phase transformations on the classical 
field configurations, just as explained after ()5. 18j) for one harmonic oscillator on phase 
space. 

We conclude that the above construction indeed yields a complex- valued wave func- 
tion <p(a), a = 1, . . . ,L, defined on the classical field configurations {-A^, . . . , A^} of 
the subsystems. It is uniquely determined up to the local phase transformations (I5.22p . 
These results make it possible to approximate a general state of the bosonic Fock 
space by our decoherent subsystems, as the following consideration shows: According 
to Lemma I5.lt the dynamics of a given bosonic Fock state can be described by a 
complex- valued wave function (p on the classical field configurations. By considering 
similar to ()5.17p a sequence of systems where the number of decoherent subsystems 
tends to infinity, we can approximate <p by a finite collection of classical field con- 
figurations {A^ 1 ' , . . . ,A( L '} and a corresponding complex-valued functions <p(a). By 
suitably adjusting the phases of the sea states ip[ a \ ■ ■ ■ ,ipf (for given reference sys- 
tems $), we can arrange that the function (p{a) satisfies (|5.2ip . Then the wave func- 
tions tpi, . . . , ipf of the whole system encode the classical potentials {A^\ . . . ,A^} 
as well as the complex- valued function (p(a), which together approximate the given 
bosonic Fock state. 

5.2. Describing a Second Quantized Fermion-Boson System. We now combine 
the considerations of the previous section with the constructions of Chapter[3]to obtain 
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a unified framework for describing second quantized fermions and bosons. We again 
consider L decoherent subsystems in the continuum limit (|5.2p . According to (|2.12p . 
we can split up the fermionic projector into the particle- and anti-particle- as well 
as the sea states. We begin for clarity in the situation without pair creation where 
the numbers n a and rif of particles and anti-particles are constant and coincide in all 
subsystems (this constraint will be removed below). Then setting re = n a — n/, the cor- 
responding many-particle wave functions of the subsystems (|4.2ip can be decomposed 
as 

= (VP A • • • A A cftf A • • • A 0g) A [Vn+i A • • • A V/] • (5.23) 

Here the round brackets can be regarded as the fermionic wave function of the particles 
and anti-particles. As explained in Section 14. 5} measurements involve superpositions 
of these many-particle wave functions, so that it is reasonable to regard the round 
brackets in (I5.23P as a vector in the fermionic Fock space 3~ eff . Likewise, the square 
brackets in (I5.23P describe the sea. The construction (I5.2ip yields a corresponding 
complex wave function <f>{a) on the classical field configurations, which can be used 
to describe the dynamics of a second quantized bosonic field. In this way, we have 
extracted from the fermionic projector both a fermionic and a bosonic quantum field. 

We now give a construction which avoids the splitting of the many-particle wave 
function into the particle/anti-particle component and the sea component. Apart from 
being simpler and cleaner, this construction has the advantage of working just as well 
for fully interacting systems, including pair creation or annihilation processes. Recall 
that in (I5.2ip we compared the many-particle wave function \E'( a ) of our subsystem 
with the wave function ^ of a "reference system" having the same classical field con- 
figuration. The proportionality factor 4>{a) then gave us the desired complex-valued 
function cf> on the classical field configurations. Giving up the requirement that the 
vector space H c iass should be represented by complex- valued functions, one can work 
instead of 0(a) with the corresponding vector 6 9yF. This has no effect on su- 
perpositions, because the complex coefficients 0(a) and 0(6) can be linearly combined 
only if the corresponding classical field configurations A^ and coincide. But then 
the corresponding Fock vectors ^W, £ 3~^ ff are linearly dependent, so that taking 
their linear combination has the same effect as taking the linear combinations of the 
coefficients 0(a) and (j)(b). This leads us to replace the complex- valued function 0(a) 

reS 

f 



constructed in (I5.2ip by a mapping with values in 3^ ff 



: {1, . . . , L} 3f : a ^ ^ A • • • A ipf . (5.24) 

In the setting involving particles and anti-particles (|5.23|) . this mapping has the nice 
property that it involves at the same time the fermionic wave functions of the particles 
and anti-particles. In free field theory, it can be thought of as the tensor product of 
a fermionic and a bosonic Fock state. As desired, two such tensor states are linearly 
dependent only if both the fermionic and bosonic parts are. Superpositions of these 
tensor states can be justified exactly as explained in Section T4.51 In a fully interacting 
system, the mapping (I5.24p can no longer be decomposed into a fermionic and a bosonic 
part, in agreement with the fact that in interacting quantum systems the bosonic and 
fermionic Fock spaces are coupled together and "mixed" by the Hamiltonian. Even in 
this in general very complicated situation, the mapping gives a conceptually simple 
description of the whole system. 
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5.3. Remarks and Outlook. To avoid confusion, we point out that the constructions 
in this chapter are not equivalent to the canonical quantization of the bosonic field. In 
particular, we do not get the physical equations for second quantized fields. Instead, we 
merely show that the dynamics of free second quantized bosonic fields can be mimicked 
by an ensemble of decoherent subsystems, each with a classical dynamics. However, 
we do not get a justification nor an explanation for the physical assumption that 
electromagnetic wave modes should behave like quantum mechanical oscillators. But 
we show that this assumption is not in conflict with the framework of the fermionic 
projector. In particular, it is possible to describe entangled bosonic states. 

In order to explain why we do not even attempt to reproduce the physical equa- 
tions for second quantized fields, we now briefly outline how interacting quantum field 
theory should be formulated in the framework of the fermionic projector. Recently, 
this formulation of quantum field theory has been worked out in detail for a system 
involving an axial field [8|. The general strategy is as follows. Instead of quantizing the 
classical field equations, we describe the interaction and the dynamics of the system 
by the action principle (j2.7f) . In the continuum limit, the Euler-Lagrange equations 
corresponding to this action principle give rise to the Dirac equation coupled to clas- 
sical bosonic field equations (cf. p. lip and (|2.13j) ) . Treating this coupled system of 
nonlinear partial differential equations in a perturbation expansion gives rise to all 
the Feynman diagrams of perturbative quantum field theory (see [U Section 8.4]). In 
particular, this gives agreement with the high-precision tests of quantum field theory. 
We remark that we get additional small corrections to the field equations which are 
absent in perturbative quantum field theory; the interested reader is referred to [SI 
Section 8.2 and 8.3]. 

Since the quantitative aspects are respected, it remains to explain the particular 
effects of quantized fields. The present paper is concerned with entanglement and 
shows that entangled fermionic and bosonic states can be described in the framework 
of the fermionic projector. For other quantum effects related to the measurement 
problem and the wave-particle duality see also [4, Section 4]. Putting these results 
together, it seems to us that the framework of the fermionic projector is in agreement 
with all effects and predictions of quantum field theory (except for the additional 
corrections discussed in [HI Sections 8.2 and 8.3]). But this agreement cannot be stated 
in terms of a mathematical equivalence, partly because standard quantum field theory 
at present has no fully convincing mathematical formulation. Also, many difficulties of 
quantum field theory clearly remain unsettled in our approach. Thus many conceptual 
and technical issues need to be debated in the future. 

We finally point out that the agreement with free quantized fields in Section 15.11 
is obtained only in the limit when the number of subsystems tends to infinity. Thus 
even for describing the quantum oscillations of the harmonic oscillator corresponding 
to one mode of the electromagnetic field, one needs to consider a large number of 
decoherent subsystems. Although this seems possible in principle, it seems hard to 
imagine that decoherence should really lead to a "fragmentation" of space-time into 
many disjoint regions with an independent dynamics. This raises the question whether 
the microscopic mixing of decoherent subsystems might not be a too simple picture for 
understanding the mechanisms of space-time on a small scale. Indeed, it might be more 
appropriate to replace this picture by a more general concept which we now explain 
in words. Recall that in Section f4.4l the decoherence of subsystems was introduced by 
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inserting a unitary transformation into the fermionic projector, 

/ 

P(x,y) = - £ U jk \4>f\x)y-<4,£\y)\ , (5.25) 
j,k=i 

if x and y are in different subsystems (see (|4.25|) ) . Since P(x,y) is a 4 x 4-matrix, a 
dimensional argument shows that there is a large class of operators U £ SU(/ ) which 
do not affect the form of P(x, y). In order to make use of this additional freedom, we 
replace U by a family of local unitary transformations U(x) 6 SU(/), which brings the 
fermionic projector to the more general form 

/ 

P(x,y) = - U ]k (x)U kl \y)\^(x)y^My)\- (5-26) 
j,k,l=i 

By dividing M into subregions and choosing U(x) to be constant on each subregion, we 
get back to the setting of Section \4Al But if U(x) is not a piecewise constant function, 
the situation is more involved. Namely, for any space-time points x and y, it is possible 
that the transformation Ujk(x) U^fy) has no effect on P(x,y) (similar as discussed 
in (|5.25p ); in this case the pair (x, y) is said to be coherent. Another possibility is that 
the transformation Ujk(x) U^ l {y) leads to cancellations in the sum so that P(x,y) 
is very small (similar as explained after (|4.25p ). in which case the pair (x,y) is said 
to be decoherent. This notion of decoherence again gives a relation between space- 
time points. But in contrast to the situation in Section 14.41 this relation is no longer 
transitive; for example, it is possible that the pairs (x, y) and (y, z) are coherent, but 
the pair (x, z) is decoherent. As a consequence, decoherence no longer gives rise to a 
decomposition of space-time into subregions. But for any fixed space-time point x, one 
can form the set M(x) of all space-time points which are coherent to x. This set can 
then be divided into subsets Mj(x) by the condition that any two points y, z £ Mj(x) 
should be coherent to each other. On the sets Mj(x), one can then again consider 
the continuum limit to obtain for example the Dirac-Maxwell system (|2.1ip . (|2. 13[) . 
Thus on the coherent space-time points one again gets a description involving classical 
field equations. Decoherent pairs of space-time points, on the other hand, are not 
connected by our action principle. We remark that it is also conceivable that two 
space-time points are partially decoherent in the sense that there are cancellations 
in the sums (|5.26p . but without P(x,y) being very small. We expect that such a 
partial decoherence would yield a large contribution to the action and should thus be 
avoided by our action principle. Then the resulting structure resembles the situation 
in Section ^. 4l in that we obtain decoherent subsystems with an independent dynamics. 
The main difference is that the subsystems are no longer localized in disjoint regions 
of space-time. Instead, they are all delocalized, and only when picking a pair of space- 
time points (x,y), the phases in the sum (|5.26p determine to which subsystem the pair 
belongs. Due to the obvious analogy to a hologram, we refer to this concept as the 
holographic superposition of subsystems (but it does not seem to be directly related to 
't Hooft's holographic principle). The main advantage of a holographic superposition 
is that a large number of subsystems no longer leads to a "fragmentation" of space- 
time into disjoint space-time regions. On the other hand, all the effects considered in 
this paper can be described just as well by decoherent space-time regions. Therefore, 
the holographic superposition is not essential for our purposes, and we shall not enter 
the detailed constructions here. 
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Appendix A. Settings where Entanglement is Impossible 

A.l. Restriction to a Subsystem. In all realistic situations, the relevant measure- 
ments are performed only in a small subsystem of the universe. Restricting attention to 
the subsystem, the effective fermionic state is no longer of Hartree-Fock type, giving 
the hope that this might make it possible to describe entanglement. The observa- 
tion that the restriction to a subsystem gives more freedom to describe the effective 
fermionic state was already used in [3j Appendix A] to show that if only one-particle 
measurements are made, describing the system by a one-particle projector P is equiv- 
alent to the general Fock space formalism (this result corresponds to Proposition IA. 51 
below). Here we generalize the methods to many-particle observables and show that 
the restriction to a subsystem does not make it possible to describe general entangled 
states. 

In order to describe the subsystem, we assume that the one-particle observables 
act only on a proper subspace 3 C IK (the "inner system"). Extending such an 
observable O by zero to all of H, we obtain a self-adjoint operator which vanishes 
on the orthogonal complement of 3. Extending this operator by (|3.13p to the Fock 
space, we obtain an operator which preserves the number of particles and vanishes on 
the orthogonal complement of the space generated by 3, 

O : 3~ n — > 3~ n self-adjoint and 0\ < A n (3«) > ± = (A-l) 

(for notational convenience the superscript of has been omitted). In order to allow 
for many-particle observables, we can generalize O to be a general operator on 3~ which 
satisfies (lA.lj) and thus leaves the number of particles n fixed. 

We again denote the projector on the corresponding /-particle Hartree-Fock state ^> 
by P f (see (pTTTjh (pT8j) and Proposition E2J) , 

Py> = /! where * = ipi A • • • A tpf . (A.2) 

Setting O = 3 1 - (the "outer system"), we decompose all one-particle states ipi into 
their inner and outer parts, 

ipi = $ + iff with ^ G 3, ip® G O . (A.3) 

Substituting this decomposition into (1A.2|) and multiplying out, one gets a sum of 
terms involving wedge products of the ijJj and ijjf ' . In order to keep track of the 
combinatorics, it is convenient to denote by I a multi- index 

I = (h,...,i g ) with 1 < ii < i 2 < ■ ■ ■ < i g < f , (A.4) 

to set \I\ := g, and to define the sign of / as the sign of a permutation in {1, ...,/} 
which maps I to the set {1, . . . , g}, i.e. 

sign(I) = (_l) i i+-+ i 9 + M ^ ii . (A.5) 
Furthermore, we denote the complement of / by O, i.e. 

= (ji,...,i h ) with 1 < ji < ••• < j h < f, g + h = f, i k ^jiMk,l. 
Finally, we set 

* / = ^ 1 A---A^ fl and *° = ^ A • • • A ^ . 
Using this notation, we obtain 

^ = sign(J) SS? 1 A , 
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and thus we can write (|A.2|) as follows, 

P f = ^2f\ sign(I) sign(l') A A ^°'\ . (A.6) 

The inner system 3 is conveniently described by a density operator, which we now 
define. A state $ G 3" (no matter if factorizable or entangled) is referred to as a pure 
state. Due to a limited knowledge on the physical system, the fermions can often not 
be described by a pure state, but merely by an ensemble of states, coming with certain 
probabilities. More precisely, one considers a family of pure states . . . , 6 3" 
together with corresponding probabilities pi, . . . ,pl, normalized as follows, 

L 

||ty||gr = l and 0<pj<l, ^2pi = l. 

i=i 

The expectation value of an observable O is then defined by 

L 

(0) = J> <*z|0tf,) . (A.7) 
i=i 

Introducing the operator 

L 

i=i 

the expectation value can be expressed in analogy to ()3.17|) by 

(O) = Tv(pO) , (A.8) 

whereas the conditions (|A.7p become 

p > and Tr(p) = 1 . (A.9) 

The operator p is referred to as the density operator. If p is a projector, it follows 
from (|A.9P that p has rank one, and thus it can be written in the form (|3.1ip with a 
pure state If p is not a projector, it is said to describe a mixed state. 

When computing expectation values of operators localized in our subsystem, we 
can take the partial trace over O to obtain an equivalent description of our quantum 
system by a density operator defined in the subsystem. This is made precise in this 
next lemma. 

Lemma A.l. For any operator O on 'J of the form (lA.lj) . the expectation value (13, 17ft 
can be expressed by (O) = Tr(pO), where p is the density operator 

f 

P = J2 E 9Kf - 3)!sign(/)sign(/') (*°'|*°) |^)(^'| . (A.10) 

3=0 I,I'with\I\=\I'\=g 

Proof. Applying (|A.1|) . we find 

(¥' A O (¥ A *°)) = A *°'| (OV 1 ) A *°) 



<?!(/-<?)! 



(*°'|*°) Tr (l* 7 )^ 7 '!©) . 



Using this relation in (|A.6P gives the result. □ 
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We point out that the density operator (|A.10|) involves states of a variable number 
of particles g = 0, . . . , /. The coefficients depend on the inner products ) of 

the wave functions in the outer region. Since in the outer region no measurements are 
possible, we cannot determine the wave functions ip®. Thus we can take the point of 
view that these wave functions can be chosen arbitrarily. At first sight, this seems to 
give a lot of freedom to choose the coefficients (\I /0 \^°) in (|A.10|) . and thus one might 
conjecture that with (|A.10j) it should be possible to describe a general entangled state. 
However, it is shown in Appendix lA.il that this conjecture is wrong, basically because 
many degrees of freedom in choosing the wave functions tp® drop out when carrying 
out the sums in (lA.lOj) . We thus conclude that the restriction to a subsystem does not 
make it possible to describe entangled states. 

Let us analyze the density operator of Lemma IA.lt In order to bring the inner 
products (^°\^° ) into a more convenient form, we introduce the (/ x /)-matrix A 
by 

Mj ■= (tfW?) ; (A.11) 

it can be considered as the Gram matrix of the one-particle wave functions in the 
outer system. This matrix positive semi-definite matrix, but this property will not be 
used in what follows. This has the advantage that our results apply even in a setting 
where the inner product (.|.) is not positive definite (for example, one might think of a 
situation where the measurement scalar product is indefinite on the microscopic scale) . 
The matrix A may be singular. Therefore, it is preferable to add a small multiple of 
the identity matrix to obtain an invertible matrix. Thus for any e > we set 

A £ = A + el and B £ = {A 6 )' 1 . (A.12) 

By Bj, j we denote the (g x g)-matrix obtained from B £ by deleting all rows and 
columns except for those corresponding to / and I', respectively. 

Theorem A. 2. The expectation value (|3.17|) can be expressed by 

(0)=Tt( p O), 

where p is the density operator 
f 

P = J2 E ff! lim^det^detCBf^))!* 1 )^'!, (A.13) 

9=0 I,I'with\I\=\I'\=g £ 
and the matrices A 6 , B £ are defined by (jATTTT) and (lA~T2j) . 
Proof. Setting h = f — g, the inner product in (|A.10|) is computed by 

<*°>°> = ^ E ® • • • ® ^l (h) I € ® • • • ® O = ^ det (A?',o) , (A.14) 

where in the last step we used (|3.ip and the definition of the determinant. Since the 
determinant is polynomial and thus continuous, it is obvious that 

det(A / ) = lim det(^Q, Q ) . 

e\0 ' 

Using this relation in (|A.14j) and (jA.lOp . we conclude that it remains to prove for any 
invertible matrix A the identity 

det(A > j0 ) = sign(J') sign(I) det(A) det(B rj ) , (A.15) 
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which relates the minors of A to the minors of its inverse. For proving (|A.15p . we 
first note that in the special case where I and I' consist of only one index, this is 
the well-known Cramer's rule for the inverse of a matrix. In the general case, the 
identity (|A.15P is stated in [131 Section 0.8.4]. It can be proved as follows. The signs 
in ()A.15P can be understood from the fact that if rows or columns of the matrices A 
and B are permuted without violating the ordering of the multi-indices I, I', O and O', 
then every such conjugation flips the sign of det(j4), and also the sign of one of the 
functions sign(J) or sign(J'). With such conjugations we can arrange that 

I = 1' = (!,.. .,g) and O = O' = (g + 1, . . . , /) (A.16) 

(but of course, the matrices A and B will no longer be Hermitian). It remains to show 
that in this case, 

det(^ ',o) = det(A) det(B rj ) . 

Using Laplace's formula, in case (|A.16|) the minor det(Ao> ,o) can be written as a 
multiple derivative of det(A), 

d d 

det(A 0/t0 ) = — detA. (A.17) 

dAu dA gg 



Using the standard formulas 

d 



dAu 



det{A) = det{A) B { 



d - B 3k = -£r^) jh = -(b (-£-A )B\ -D P D 



ik 

jk 



dA u J ^ dAa y V \dAii 

one can iteratively carry out the derivatives in (|A.17|) to obtain 
det(Ao',o) = det(.A) x (homogeneous polynomial of degree g in Bjk with j, k € J). 



Going through the combinatorial details, one finds that this homogeneous polynomial 
coincides precisely with det(Bpj). □ 

We now illustrate this theorem by a few examples and work out simple consequences. 
We first explain how to recover a Hartree-Fock state. 

Example A. 3. (Hartree-Fock state) Choosing vp** = ^ and = 0, the Gram 
matrix (|A.11|) vanishes and thus 

A £ = e and B e = e' 1 . 

Hence 

det(A e )det( J B|, / ) = e f ~^ 5 I>P , 
and this vanishes in the limit e \ unless I = (1, . . . , /). Hence (|A. 13|) reduces to 

P = /! !#)(*!, 

so that the density operator coincides with the projector Pf ()A.2|) of the whole system. 



Example A. 4. (A mixed state) We choose / = 3, ipf = 0, = \\i/j§\\ = k 6 [0, 1] 

and = 0. Thus 

A e = diag(e, k 2 + e, k 2 + e) , 
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and a short calculation yields that 

p =k>?>(V?I + 2k 2 A ^1)^1 A $\ + |V? A A ^|) 

+ 3!|^AV 2 3 AV 3 3 )(^A^A^|. 

Thus the density operator of the subsystem is an ensemble of one, two and three 
particle states. The one-particle state V>i is always occupied; this is because the cor- 
responding eigenvalue of the matrix A vanishes. The two-particle component cannot 
be represented by a pure state; it is a mixed state. 

We next show that when restricting attention to one-particle observables, the ex- 
pectation value of any vector in the fermionic Fock space can be approximated to 
arbitrary precision by a one-particle projector P. 

Proposition A. 5. For every normalized n-particle state \P in our subsystem, i.e. 

* € <A n (3 n )> C J n and ||*|| = 1 , 

there is a sequence of projectors (-Pfc)fceN * n ^ of finite rank such that for every 
one-particle observable O, 

(tf = lim Tr(P k O) . 

k—>oo 

Proof. Writing $ as a linear combination of Hartree-Fock states and computing the 
expectation value with (|3. 13[) and (|3.ip . we obtain 

k,l 

with suitable vectors € 3 and complex coefficients c k i- Using an approximation 
argument, it clearly suffices to consider finite sums and finite linear combinations. 
Choosing an orthonormal basis ipk of the subspace (ipj*,) C 3, and expressing the as 
linear combinations of the t/jf-, we obtain 

k,l 

Diagonalizing the symmetric matrix pi-i by a unitary transformation, this representa- 
tion simplifies to 

=Y,PkWk\Oil> k ). (A.18) 

k 

Since the operator is positive and normalized, it follows that 

p k > and y Pk = l. 

k 

We choose an orthonormal family of vectors (4>k) in O and set 

Vfe 0t = \/~Pk^k + V 1 - Pk 4>k ■ 
Then the family (ip^) is orthonormal, and thus 

p=^>n<vri 

k 

defines a projector in "K. A short calculation shows that the expectation value Tr(PO) 
coincides with the right side of (|A.18p . □ 
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The next example shows that the restriction to a subsystem does not make it possible 
to describe general entangled states. 

Example A. 6. (The spatially separated singlet state) Let us try to realize the spatially 
separated singlet state of Example 13.11 We assume that the inner system 3 C "K has 
the orthonormal basis (ip\, ipg, ip^, V"b)- Our goal is to find a projector P in J£ such 
that the corresponding density operator of the subsystem (|A.13|) coincides with the 
projector onto the singlet state (|3.3f) . i.e. 



\i>l A 4xv>i a 4i - 14 a a 4 



(A.19) 

The simplest way to see that such a projector P does not exist is to observe that 
in (|A.13P one necessarily gets the contribution involving four particles 



\if>\ A il>% A ifa A A V B A A ^| , 

which is not present in (1A.19P - However, this argument is not fully convincing, because 
an additional four-particle contribution would not be observable in the standard spin 
correlation experiments where Alice and Bob can only detect one particle at a time. 
Furthermore, by making the matrix elements of B £ small, one could try to arrange 
that the four-particle contribution is negligible. For these reasons, it is preferable to 
show that we cannot even realize that the two-particle component of (1A.13P coincides 
with ([AT9]) . as we now explain. 

Restricting attention to the two-particle component of (IA.13j) . we need to consider 
the determinants of the 2-submatrices of the matrix B £ . Noting that in the case |/| = 
1 1'\ = 2, we can write the product of the determinants in ()A.13P as 



det(yl £ )det( J B|, 7 ) = det ^det{A £ ) B e V I 
it suffices to consider the 2-minors of the matrix 



B 



B := lim y/deb(A e ) B e . 

Grouping the first two and the last two basis vectors together, we write B as the block 
matrix 

Bn B\2 

Bi\ #22, 

whose entries are 2x2 matrices acting on the subspaces <{ip\, ^s}> an d <V^> yb^ > 
respectively. In order to realize (|A.19p . we need to arrange that the determinants 
of the submatrices fin and B22 are non-zero, but all other 2-minors must vanish. 
Diagonalizing B\\ and B22 by unitary transformations in the subspaces <{tp\,ip B }> 
and <{ip^,ip B }> , respectively, we obtain 



B 



/pi a c\ 

p 2 b d 

a b P3 

y c d pi j 



with pi > 



Evaluating the condition det Bjji = for / = (1,3) and I' = (1,2), we see that b 
must vanish. Similarly, taking / = (1,4) and I' = (1,2) yields d = 0. Repeating 
this procedure for / = (2, 3) and / = (2, 4) we find that a = c = 0. Hence B is a 
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positive diagonal matrix. But then the submatrix Bj j with / = (1,3) has a non-zero 
determinant, a contradiction. 

Remark A. 7. (Systems with purely classical bosonic fields) Let us consider a system 
where second-quantized fermions are coupled to classical bosonic fields, in the pure 
sense that also all measurement devices can be described by the classical field equations. 
As pointed out in [3, Appendix A], the coupling of the fermions to classical fields can be 
described by one-particle observables (like the expectation value of the Dirac current 
in Maxwell's equations or the expectation value of the energy-momentum tensor in 
Einstein's equations) . In view of Proposition lA.5l the resulting system can be described 
by a one-particle projector P. On the other hand, we saw in Example IA.6I that this 
framework does not allow for the description of general entangled states. We conclude 
that the physical phenomenon of entanglement makes it necessary to consider second 
quantized bosonic fields. 

This consideration suggests that there is a close connection between entanglement 
and second quantization of the bosonic fields. This connection becomes clearer in 
Chapters H] and O because the method of microscopic mixing, which allows for the 
description of entanglement (see Chapter H]), also makes it possible to describe second 
quantized bosonic fields (see Chapter [5]). 

A. 2. Microscopic Mixing of the Wave Functions. In order to analyze the setting 
of Section 14.11 in more detail, it is useful to get a connection to the fermionic Fock 
space formalism. For a one-particle wave function dp 6 Ji, the creation operator a^(da) 
is defined on Hartree-Fock states by 

a) {do) (V>i A • • • A if> n ) = do A Vi A • • • A Vn • 

By linearity, it is extended to an operator a\do) : 3~ ra — > S'n+i- Its adjoint a(do) : 
9~n+i — > 3>i> the so-called annihilation operator, acts on Hartree-Fock states by 

n+l 

a(4>) (ipi A ■ ■ ■ A Vn+i) = ^( _1 ) fc+1 ^1^) ^1 A • • • A Vfc-i A Vfc+i A • • • A t/i n+1 • 

k=l 

A straightforward calculation shows that the creation and annihilation operators sat- 
isfy the anti-commutation relations 

{a(0),a^)} = (#/>>!?■ (A.20) 

A one-particle observable (|3.13p can be expressed in terms of the creation and annihi- 
lation operators by 

k,l 

where {dad) denotes an orthonormal basis of (!K, (.|.)). Products of one-particle observ- 
ables can be transformed with the anti-commutation rule (|A.20|) ; for example, 

OlO%= » t (^fci)a t (^fa)^fci|Oi0i x )(0 fa |O 2 ^ a >a(^ l )o(0, a ) (A.21) 

ki,k2,h,h 

+ J2<f(<j> k ){<l> k \0 1 2 <l n )a(<fn), (A.22) 

k,l 

where in the last line we used the completeness of the basis (4>i). 
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A useful rule in quantum field theory is Wick ordering, denoted by colons, which 
states that all creation operators should be written to the left and all annihilation 
operators to the right, leaving out all terms which would be generated by the anti- 
commutations. For example, Wick ordering the above product of one-particle observ- 
ables amounts to omitting the term (|A,22|) . 

:OfOl-.= Yl ^(^J^^J^lOl^)^!^^}"^)"^)- (A.23) 

A general Wick-ordered two-particle observable can be written as 

:0:= Yl a) {(j) kl ) a) {(t) k2 ) g{h, k 2 , h,h) a(^)o(fe), (A.24) 

where the function g is anti-symmetric in its first and last two arguments, 

g(k 1 ,k 2 ,hj2) = -g(k2,h,h,h) = -fl(fci, h, h, h) ■ (A.25) 

In the next proposition we express the expectation value of this two-particle observable 
in terms of the one-particle projector P. 

Proposition A. 8. Describing a many-fermion state by a projector P in (JC, (•!•)), 
the expectation value of the two-particle operator (|A,24p with g according to (|A.25|) is 
given by 

(:0:) = Y (i&k\P4>k){<i>i\P<h) ~ mP<Pi){4>i\P^k)) g(k,l,k,l) . 

kytl 

The expectation value of the Wick ordered product (|A.23P is 

(:Of0f :) = Tr M (PC 2 ) - Tr % (PO l P0 2 ) . (A.26) 

Proof. We again represent the projector Pj in the form (|3.1ip . Using the obvious 
transformation law of the function g under basis transformations, we can arrange 
that 0i = tpi, . . . , <pf = ipf. Then the results follows from a straightforward calcula- 
tion. □ 

Note that the formula (|A.26|) coincides with (|3.15p . except that the summand Tr^PCiC^) 
is now missing as a consequence of the Wick ordering. 



Proposition A. 9. In the setting of Section \4-l\ there is no projector P in Ji which 
reproduces the expectation values of the spatially separated singlet state (|3.3p with re- 
spect to the spin operators of Alice and Bob and the corresponding two-particle spin 
correlation operators. 

Proof. We let 5| be the spin operator having the expectation values one if the spin 
is up and zero if the spin is down. Similarly, S± is the operator for spin down. The 
spin operators corresponding to the observers Alice and Bob are denoted by S-ma 
and Sy^B, respectively; they are operators in "K. Taking for convenience the Wick 
ordered products of the corresponding operators on the Fock as defined by (|A.23p . we 
can compute the following expectation values of the singlet state (|3.3|) , 

(Sf,A) = (Si,b) = 
1 2 (A.27) 

(: Sy A Sf^ B :) = - , (: S^ A S^ A :) = (: Sf B Sf B :) = 
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(all the other expectation values are irrelevant for the proof). Assuming that the 
statement of the proposition is false, these expectation values can be reproduced by a 
suitable one-particle projector P. Applying Lemma 13.31 and Proposition IA.8( we can 
then express the expectation values by traces over the one-particle Hilbert space "K, 

(0?) = Ti(PO) and (:OfOf :) = Tt{PO x ) Tr(P0 2 ) - Tr(P0iP£> 2 ) . 

Using the idempotence of P, the arguments of these traces can all be rewritten purely 
in terms of the operators 

T %A :=PS^ A P and T UB := PS itB P , (A.28) 

which are symmetric and of finite rank. The relations (IA.27P give rise to the conditions 

Tr(T t>A ) = Tr(T i>B ) = \ , Tr(T t 2 A ) = Tr(T^ B ) = I , Tr(T t>A T i>B ) = ~ . (A.29) 

At this point it is helpful to regard the operators T . as vectors in the real Hilbert space 
of symmetric Hilbert-Schmidt operators with the scalar product (A,B)hs = Tr(AB) 
and corresponding norm ||.A[|hs = \J (A, A)hs- Then the last two equations in (|A.29|) 
imply that 

\\T t , A + T i>B f HS = Tr(T 2 A ) + 2 Tr(T tiA T 4jB ) + Tr(T^ B ) = . (A.30) 

It follows that T^ )A = —T^Bi m contradiction to the first equation in (|A.29p . □ 

We point out that in this proof we did not use that the one-particle operators are 
invariant on the subsystems (|4.4p . Thus dropping this assumption would not change 
the statement of Proposition IA.91 

Acknowledgments: I would like to thank Andreas Grotz, Mario Kieburg, Joel Smoller, 
Roderich Tumulka, Tilo Wettig and the referees for valuable comments on the manu- 
script. 
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